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Using the lattice QCD equation of state for an isentropically expanding fireball we follow the 
evolution of the effective number of particles in an ideal gas, N c g — pV/T. This number reduces to 
its third around the crossover temperature, which helps to resolve the entropy problem inherent in 
hadronization models assuming quark coalescence. 



The hadronization of quark matter is most probably a complex, highly non-equilibrium process. During a sudden 
cooling and expansion elementary particle correlations are re-arranged while passing the crossover temperature of 
QCD. It is therefore surprising that simple models, such as a nearly ideal quasiparticle gas in quark matter and a 
^NJ ■ hadronic resonance gas, are good fits to the equation of state of the strongly interacting system studied in lattice QCD. 

Moreover during a sudden, quark-recombining hadronization scenario the number of effective degrees of freedom is 
Qh dramatically reduced. 

Hadronization models therefore cither dynamically fragment quarks and gluons to several hadrons, or face the 
problem how not to reduce the entropy. The simple quark-coalescence idea, based on the valence quark model of 
hadrons, redistributes and bounds quasiparticles into hadrons assuming an effectively low gluonic quasiparticle number 
in the formation process. The question naturally arises whether a small quasiparticle number can be realistically 
assumed while not reducing the entropy of the quark matter. The argumentation pointing out an entropy reduction 
usually assumes a massless ideal gas full with gluons on the QGP side and only the lightest hadrons (pions) on 
the hadron side. This is not very realistic. We show that the entropy obstacle can be circumvented by choosing 
appropriate temperatures on the quark and hadron matter sides in accord with the equation of state obtained from 
; lattice QCD. 

In this paper, supported by our recent analysis of lattice QCD equation of state (eos) data[l|, we show that during 
an adiabatic expansion the number of massive particles in a quasiparticle approach Q can be reduced to half or to 
one third in a realistic expansion time. This quantitative property of the lattice QCD eos explains the success of the 
massive quark matter coalescence hadronization model[j|, 0, Q. 

The following fit to the continuum extrapolated pressure curve of Ref . [fj can be considered (cf . FigQ]) 

p = cT i a(T c /T) (1) 
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°(z) = ^—T — nre~ Az , (2) 



where z = T c /T,a = 0.91,6 = 0.11 and A = 1.08. The Stefan-Boltzmann coefficient per degree of freedom, cr(z), is 
temperature dependent. This reflects the interaction, in particular mass generation, in the lattice QCD eos[l|. At 
infinite temperature the massless ideal gas is expected, therefore cr(0) = 1. Of course the data points of the QCD 
lattice eos also can be used without any interpolating fit. 

The number of particles in an ideal gas, which has the same pressure as a given interacting system, can be defined 
through the thermodynamical potential. We refer to this quantity as the effective number; it represents the number 
of non- or weakly interacting quasiparticles: 

N cS = ^n(V,T,[i). (3) 

For homogeneous systems, as it is assumed by the extrapolation of lattice QCD data to infinite volume, this potential 
is simply related to the pressure: f2(V,T, /i) = Fp(T,/i). The elementary particle number is on the other hand given 
by 
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FIG. 1: The normalized pressure as a function of the temperature fitted to lattice QCD eos data of Ref.[i 



It is straightforward to see that the expression in the bracket vanishes and hence N e s = N in the Boltzmann 
approximation 

Sl(V,T,fi) = n (V,T)e^ T . (5) 

In this case, using homogeneity, the Boyle-Mariotte law is fulfilled: 

P V = NT. (6) 

For the quark coalescence model of hadronization the quasiparticle number, N e g is relevant. For a mixture of 
quasiparticle gases 

f7 = ^a0/T)e c ^/ T , (7) 



and the effective number becomes 
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It is a weighted sum of elementary particle numbers. In order to estimate the hadron numbers by the coalescence 
principle, this quantity has to drop suitably while the temperature drops and the volume expands. The main question 
to clarify is whether keeping the entropy constant during the expansion an appreciable reduction in the quasiparticle 
number N e g is possible in accordance to the equation of state of an interacting quark-gluon plasma. 

Using the definition of entropy, S — dfl/dT, the effective number of equivalent non- interacting particles (quasi 
quarks and gluons above T c , hadrons below T c ) divided by the total entropy can be extracted from the p(T) curve as 
follows 
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Substituting our fit for a{z), eq.©, we arrive at 
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AT eg = 5f4 + ^ + 6r(i + e ; /& _ Te/feT) j . (10) 
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FIG. 2: The effective number of particles divided by the total entropy for a non-interacting ideal gas, which reconstructs lattice 
QCD pressure as a function of temperature. A reduction to 1/3 is achieved by cooling from T2 = 254 MeV to Ti = 137 MeV. 

Figd] plots the ratio N c g/S as a function of the scaled temperature T/T c using the parameters a = 0.91,6 = 
0.11, A = 1.08. The big dots indicate the result obtained from p(T) data by numerical derivation. 

The resulting curve in Fig [2] shows two turning points: one below, one above T c . Considering a standard value, 
T c = 167 MeV, the turnings are located at Ti = 0.82T C = 137 MeV and T 2 = 1.52T C = 254 MeV. These points were 
estimated as the endpoints of the steepest linear rise in the N c s/S curve. The effective number of ideal quasiparticles 
falls to its one third between these two temperatures while keeping the total entropy constant. This roughly corre- 
sponds to changing from a colored, quark dominated gas to a color singlet one. The temperature T 2 is around the 
color deconfinement temperature obtained earlier in pure gluon systems [7|. 

The fastest change in N g/S is concentrated to the temperature range between T\ and T 2 . This may be considered 
as the transition region between the high temperature QGP (where the quasiparticle masses are set by the scale 
XT C = 180 MeV) and the hadronic phase near T c (dominated by heavy resonances [8j). The large number of different 
hadron types accommodates the high entropy inherent in the QGP with liberated color degrees of freedom. 

In an earlier work[l| we have analyzed the lattice QCD equation of state in terms of distributed mass ideal parton 
gas. We found that at temperatures over T c the mass distribution starts at XT C w 180 MeV and has its maximum 
around the valence quark mass, m v « 300 — 350 MeV. Below T c the characteristic minimal mass is about ten times 
higher, 10T C w 1.67 GeV. 

We interpret Fig[2] in the following way: as the temperature decreases due to physical expansion of the volume at 
constant entropy, the number of quasiparticles drops by about a factor of three. We know from earlier studies (l3| 
that around the same temperature the characteristic mass scale of the distributed mass quasiparticles can increase by 
a factor of ten. This reduction in the effective particle number definitely does not contradict to the constant entropy. 
Thus we conclude that hadronization models utilizing such a reduction in the number of particles do not necessarily 
reduce the total entropy. 

In order to connect this observation to the real time behavior of expanding quark matter one has to explore 
the relation between the eos and the hydrodynamical expansion of the perfect fluid. This expansion is isenthalpic, 
dE + pdV = TdS + ^i/iidNi = 0, if chemical transmutation rates are neglected even isentropic dS — 0. Since 
the QCD color deconfinement transition is a crossover, no latent heat is released, and therefore the hadronization is 
neither an exothermic nor an endothermic process. We consider dS — 0. A local volume is defined by the expansion 
rate 




(11) 



where d/dr = u^d^ is the comoving derivative. The cooling takes place due to the reduction of the internal energy 
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FIG. 3: The adiabatic cooling of QCD matter as a function of the expanding volume. The horizontal lines belong to the 
temperatures T2 and Ti, respectively. 



E = eV by the work of the pressure p on the local expansion: 

de , . 1 dV „ . „, 

d-r + ( e+P) V^ = °- (12) 

For a one-dimensional cylindrical expansion, V — ttR 2 t with the proportionality constant dV/dr — 7ri? 2 , the above 
equation describes Bjorken's flow: 



de 



0. (13) 



If a transverse flow is present, the change of the local volume defined by the local four-divergence of the relativistic 
velocity field (eq. [TTj) may lead to a more complicated V(r) relation. This effect can only be explored by solving 
Euler's equation, too. Some special analytical solutions without comoving pressure gradient, including transverse flow 
are given in 0. 

In our further analysis we display the cooling, i.e. the evolution of the temperature T as a function of the expanding 
volume V. This way we do not specify the time evolution of the volume factor, or the four-divergence of the flow; 
this presentation is independent of the details of the transverse flow. 

By using lattice QCD equation of state a first integral of the cooling-expansion law can be obtained by keeping 
S = Vs = Vdp/dT constant. In Figj3]the temperature is plotted against the expanding volume of a matter satisfying 
an equation of state fitted to lattice QCD simulations. The starting point was T — 2.5T C , the endpoint T = 0.8T C . The 
change from a fast cooling to a slower one occurs near T c ; the double- logarithmic plot shows a kick in the power-law 
behavior clearly. 

In FigH]the ideal particle number N e g = pV/T is plotted as a function of the expanding volume V. Both axes are 
scaled with the respective initial values. The steepest drop can be observed around the T — T c temperature region. 
It is realistic that N e g is reduced by almost a factor of three between T2 and T\ , while the volume expands to several 
ten times of its original size. We conclude from this observation that quark-coalescence models can be successful. 
The number of colors is reduced out from the effective number of degrees of freedom due to confinement even during 
a mild expansion, like an ideal fluid. Further non-equilibrium effects are likely to make the expansion faster in real 
experiments, but nowadays ab initio QCD calculations cannot yet deal with such a situation. 

In conclusion we observe that the lattice QCD equation of state leads to a reduction of about one third in the number 
of equivalent quasiparticles between the endpoint temperatures of the linear part in the N e g/S — T curve. The upper 
temperature T2 is around the phase transition temperature of pure gauge field systems, and the lower temperature T\ 
is near to the hadrochemical decoupling temperature observed in relativistic heavy-ion collisions. This may signalize 
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FIG. 4: The ideal number of particles N c f[ = pV/T as a function of the expanding volume by adiabatic cooling of the QCD 
matter. The temperature labels T2 and Ti indicate the corresponding volumes. 



that the hadronization of quark matter is dominated by recombination processes of massive quarks, as it is assumed 
in various coalescence models [3I HoL 11, 12 1. This reduction can take place at constant total entropy, contrary to some 
claims raised at diverse occasions. 
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